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ABSTRACT 



Context. In the context of exoplanet detection, the performance of coronagraphs is limited by wavefront errors. 
Aims. To remove efficiently the effects of these aberrations using a deformable mirror, the aberrations themselves must be measured 
in the science image to extremely high accuracy. 

Methods. The Self-Coherent Camera which is based on the principle of light incoherence between star and its environment can esti- 
mate these wavefront errors. This estimation is derived directly from the encoded speckles in the science image, avoiding differential 
errors due to beam separation and non common optics. 

Results. Earth-like planet detection is modeled by numerical simulations with realistic assumptions for a space telescope. 
Conclusions. The Self-Coherent Camera is an attractive technique for future space telescopes. It is also one of the techniques under 
investigation for the E-ELT planet finder the so-called EPICS. 

Key words, instrumentation: adaptive optics — instrumentation: high angular resolution — instrumentation: interferometers — 
techniques: high angular resolution — techniques: image processing 



1. Introduction 



Telescope 



Deformable mirror 



Very high contrast imaging is mandatory for the direct detec- 
tion of exoplanets, which are typically a factor of between 10 7 
and 10 10 fainter than their host and often located within a frac- 
tion of an arcsecond of their star. First of all, coronagraphs are 
required to suppress the overwhelming flux of the star but they 
are limited by wavefront errors in the upstream beam, which 
creates residual speckles in the science image. Adaptive optics 
must be used to correct for the effect of most of these aberra- 
tions. Some remain uncorrected generating quasi-static residual 
speckles ((Cavarroc et all |2006[) . Interferential techniques take 
advantage of the incoherence between companion and stellar 
lights to measur e these wavefront errors in the science image 
to high accuracy (Codona & Angel, 2004; Guy on, 2004). In this 
Letter, we describe such a technique called a Self-Coherent 
Camera (JBaudoz et all 12006). Residual speckles in the science 
image, also called interferential image hereafter, are spatially 
encoded by fringes so that we can derive an estimation of wave- 
front errors to be corrected by a Deformable Mirror (DM). Since 
the number of DM actuators is finite, this correction leaves resid- 
ual speckles. Thus, after reaching the DM limit correction, we 
apply an image post-processing algorithm (iGalicher & Baudozi 
2007) to achieve Earth-like planet imaging. Hereafter, we detail 
the SCC principle. Then, we describe the wavefront error esti- 
mator that we use. Finally, we present expected performances 
from space. 



Beamsplitter 



Image channel 
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Reference channel 
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Fig. 1. Self-Coherent Camera principle schematics. 



2. Principle and aberration estimator 

The beam from the telescope is reflected onto the DM and is 
split into two beams (Fig. [J}. The image channel (shown in red 
in electronic edition) propagates through a coronagraph. It con- 
tains companion light and residual stellar light due to wavefront 
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errors. Its complex amplitude is ^siO + ^c(^)> where £ is the 
pupil coordinate, ^s and ^c represent the stellar and compan- 
ion complex amplitudes of the field in the pupil plane respec- 
tively, just after the D diameter Lyot stop. The second beam, 
called the reference channel, is filtered spatially in a focal plane 
using a pinhole of radius smaller than A/D. Almost all compan- 
ion light is stopped since it is not centered on the pinhole. In the 
pupil plane just after the diaphragm (£>r), the reference com- 
plex amplitude is called ^r©. The pinhole reduces the impact 
of wavefront errors on ^r since it acts as a spatial frequency 
filter. An optic recombines the two channels, separated by £o in 
the pupil plane, and creates a Fizeau fringed pattern in the focal 
plane. Residual speckles are therefore spatially encoded unlike 
companions. The mean intensity of residual speckles of the im- 
age channel is almost spatially flat and attenuated by the corona- 
graph. To optimize the fringe contrast, we have to match the in- 
tensity distributions and fluxes of image and reference channels. 
We use a £>r < D diameter diaphragm to obtain an almost flat 
reference intensity in the focal plane. This diaphragm reduces 
the impact of aberrations on ^r, since only a fraction of the 
diffraction peak of size A/Dr is detected in the image (image b 
in Fig. O. This implies that the reference channel is quite insen- 
sitive to ab errations and can be calibra ted before the interference 
recording (Galicher & B audozll2007l) . Fluxes are equalized us- 
ing a variable neutral density in the reference channel before the 
pinhole (Sect.[3]). In Fig. [2 we present, on the same spatial scale, 

(a) the image formed after the sole coronagraph for a pupil of 
diameter D (the sole image channel) showing residual speckles, 

(b) the image corresponding to the sole reference channel for a 
pupil of diameter D R and (c) the interferential image, where the 
residual speckles are spatially encoded by fringes. 




Sole Image channel 



Sole Reference channel 




Science image 

Fig. 2. (a) Image formed after the sole coronagraph for a pupil 
of diameter D (sole image channel) showing the residual speck- 
les, (b) Image of the sole reference channel for a pupil of di- 
ameter Dr. (c) Interferential image (science image) where the 
speckles are spatially encoded by fringes. The spatial scale is 
the same for all the images. 



In polychromatic light, the intensity 1(a) of the interferential 
image on the detector is 



^) = Jtir[/s(f) + ik(f) + /c(f) 

+ 2^(A s (f)A*(f)exp(^))W, 



(1) 



where a is the angular coordinate in the science image, A; 
the Fourier tranform of the corresponding ¥,-, It the inten- 
sity |A/| 2 , and A* the conjugate of A/. The wavelength A belongs 
to <R = \Ap - A/1/ 2, Ao + A/l/2]. Following the work by 
Bor de & Traubl (|2006), we estimate wavefront errors from resid- 
ual speckles in the science image. For this purpose, we propose 
to extract the modulated part of /, which contains a linear com- 
bination of A s and A R . First, we apply a Fourier Transform on / 
and isolate one of the lateral correlation peaks. We then apply an 
inverse Fourier Transform and obtain /_ 



2ina£o 



dA 



(2) 



In Eq.|2j As and A£ depend on aD/A, inducing the speckle dis- 
persion with wavelength. Fizeau interfringe A/ go is proportional 
to wavelength in the exponential term. Both effects degrade 
the wavefront estimation from 7_ when the useful bandwidth 
is large. However, the fringe wavelength dependence is domi- 
nant. It may be more appropriate to consider an Integral Field 
Spectrometer at modest resolution (R = 100), or to simulate 
the use of a short bandpass fi lter with a c hromatic compensator. 
Such a device, proposed by Wynne (1979), almost correct for the 
two chromatic effects over a wide spectral band (A A ^ 0.2 Aq) 
to provide a smaller effective bandwidth (AAqs =* 0.01 Aq). It 
enables us to assume as close as possible a monochromatic 
case in our model of SCC image formation. We firstly assume 
A/leff <^ Ao, so that As and A^ are constant over the spectral 
band. We obtain from Eq.[2] 



™>"*ffl*m*-> 



linage 



dA 



(3) 



We have to estimate As or, more precisely, its inverse Fourier 
transform Ts and we deduce 



Vs(& - r~ 



I-(a)F*(a) 



(4) 



(5) 



[A* R (aD/A )\\Ff\ 
where < F~ 1 denotes the inverse Fourier transform, 

F = I — exp (2 in a go/A)dA 

and F* its conjugate. 

As a second assumption, we consider that wavefront errors 
(p we are attempting to measure are small and we can write the 
star field ^F's in the pupil plane upstream from the coronagraph 
as 



¥' s (£)^oP(£) 1 + 



2 I 7T 0(g) 
^0 



(6) 



where ¥0 is the amplitude of the star assumed to be uniform 
over P, which is the unitary flat pupil of diameter D. 

In a third step, we ass ume a perfectly achromatic corona- 
graph ( Cavarroc et all 120061) . which allows us to remove the co- 
herent part of the energy ¥0 P to ^'s 



Aq 



(7) 
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Finally, Eq.|4] and [7] provide an estimator of the wavefront errors 
within the pupil 



0(f) 



271 



I{T- 



I-(a)F\a) 



*¥ A* R (aD/A )\\F\f 



(8) 



with I{ } the imaginary part. In Eq.[8l F depends only on known 
physical parameters, £o and the spectral bandwidth, and is nu- 
merically evaluated. We can estimate ¥0 since we can calibrate 
the incoming flux collected by the telescope. T_ is derived from 
the recorded image I. Finally, we have to divide by the complex 
amplitude A R , previously calibrated. Setting D R « D, we ob- 
tain an almost flat reference intensity and therefore avoid values 
close to zero in the numeric division. We notice in Eq.[8]the lin- 
ear dependence of on /_, which is measured directly from the 
interferential image. We attempt to correct for these wavefront 
errors, estimated from Eq. [8] using the DM. Then, we record 
a new interferential image in which quasi- static residual speck- 
les have been suppressed and companions are now detectable. 
Practically, few iterations are required to reach high contrast un- 
der our assumptions used to derive the estimator and because of 
noise. 

We note that we require SCC sampling sufficient to de- 
tect the fringes that encode the residual speckles. The sampling 
is then larger than the classical s amp ling used in earlier tech - 
niques proposed bv lGuvonl (12004 and lCodona & Angell (2004). 
However, the SCC needs a single image to estimate wavefront 
errors, whereas the other two techniques require at least two im- 
ages since they use an on-axis recombining as in a Michelson 
scheme and either temporal or spatial phase shifting arrange- 
ments. Finally, instead of spreading the incoming light into sev- 
eral images, the SCC spreads the light into fringes. 

3. Performances 

We consider an SCC device operating in visible light (Aq = 
O.Sfim, AA =* 0.2 Ao, AA eS = 0.01 Ao, section©. We as- 
sume a perfectly achromatic coronagraph. The beamsplitter in- 
jects 99% of the incoming energy into the image channel. 
The filtering pinhole radius is Aq/D and D equals 25T)r. To 
be more realistic, we assume a calibration of the reference 
channel with a non- aberrated incoming wavefront and enter 
this value into the expression for the estimator in Eq. [8] We 
consider a 32 x 32 DM. The nth-actuator influence function 
is exp (- 1 .22 (32 (£ - £ n )/T)) 2 ), where £ n is the center of the nth- 
actuator. We call *H the (32 Aq/D) 2 corrected area which is cen- 
tered on-axis. We chose £0 = 1.05(1.5/) + 0.5 T) R ) to ensure 
that the correlation peaks of T(I) did not overlap, which corre- 
sponds to about 1.5 interfringes per A/D. We use 1024 x 1024 
pixel interferential images with 4 pixels for the smallest inter- 
fringe over H (Shannon criteria). Compared to the classical sam- 
pling used in Guyon's and Codona's devices of 2 pixels per A/D, 
the SCC image is 6/2 = 3 times oversampled, which reduces the 
field of view a priori. However, if the read out noise is not a lim- 
itation, this oversampling is not a problem since the interesting 
area <H is given by the number of actuators of the DM. We con- 
sider static aberrations in the instrument upstream of the coron- 
agraph. We adopt a 20nmrms amplitude with a spectral power 
density varying as f~ 3 , where / is the spatial f requency, which 
corre sponds to typical VLT optic aberrations dBorde & Traubl 
2006). We simulate an 8 m-diameter space telescope with a 50% 
throughput pointing a G2 star at lOparsec. The quantum effi- 
ciency of the detector is 50%. We consider photon noise, set the 
read out noise to 5e- per pixel, and consider the zodiacal light to 



be a uniform background at 22.5 mag.arcsec -2 . We have not used 
any linear approximation to simulate the focal plane images. 

By correcting the wavefront errors, we improve the corono- 
graphic rejection and the reference intensity T R becomes dom- 
inant in the science image / (Eq. [I]). We adjust the calibrated 
neutral density in the reference channel at each step. To deter- 
mine the value of the neutral density, we estimate the ratio r of 
the incoming energies from the image and reference channels in 
the center of the image to be r = J (1(a) - Tr(oO)/ JL, feW 
where W represents the (22Ao/D) 2 centered on-axis area, and 
optimize fringe contrast for the next step. The neutral density 
transmission is 1 at step and 5.6 x 10" 3 at step 3. At 5 A/D, 
the average number of photons per pixel is about 1.3 and 50 for 
the reference channel, and 150 and 115 for the image one, re- 
spectively at steps and 3. Finally, the intensity in *K decreases 
as the coronagraphic rejection increases. At each step, we ad- 
just the exposure time to optimize the signal-to-noise ratio in the 
16-bit dynamic range of the detector. 
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SCC post -processed image with 20nm rms (3 hours) 
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Iteration 1, Exposure time = 300s 

lteration2, Exposure time = 637s 
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Fig. 3. 5<x detection limit vs angular separation. 



We define the 5 cr detection d$ a to be 



30 



d5cr(p) = 



5(T(p) 

/o ' 



(9) 



where cr(p) is the azimuthal standard deviation of the considered 
image at the radial separation of p and To the maximum inten- 
sity of the central star without a coronagraph. We plot ds a for 
the interferential image T versus angular separation for several 
iterations of the correction (Fig. [3]) . The 5<x dete ction limit cor- 
responds to an azimuthal average (JCavarroc et all 12006). In the 
figure we specify the exposure time of each step. In iteration 0, 
we measure the coronagraphic residue due to the 20nmrms 
static aberrations without any correction. The algorithm con- 
verges in a few steps (~ 3). The dashed green line represents 
the coronagraphic image, without SCC, computed with a full 
correction by the 32 x 32 DM. This curve is almost surperim- 
posed on the curve of iteration 3. This illustrates that the SCC 
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is limited by the aberrations linked to the DM uncorrectable 
high-order frequencies. The level of this limit depends only on 
the num ber of actuators of the DM and the initial aberration 
level dBorde & Traubl 120061) . To improve the performance, we 
may increase the number of DM actuators or use higher qual- 
ity optics. In a second step, we apply to the final iteration im- 
age the post-processing algorithm that we presented in previous 
papers (Baudoz et al., 2006; Galicher & Baudoz, 2007). The 5cr 
detection limit of the SCC post-processed image is plotted in 
Fig. [3] (full red line). The increase in the faintness corresponding 
to the 5 cr detection is about 10 5 at 5Aq/D in a few steps. An 
Earth-like planet, 2 x 10" 10 fainter than its host star, is detected 
at the 5 cr confidence level in about 3 hours. Contrast outside <H 
is improved slightly during the first steps because both the refer- 
ence flux (neutral density) and the corresponding noise decrease. 
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Fig. 4. SCC post-processed image of the final iteration image 
corresponding to a total exposure time of ~ 3 hours. 2 x 10" 10 
companions are present at 1, 3, 5, and 1Aq/D on a spiral The 
field of view is about 32 x 32 A/D. The intensity scale is linear. 



achieved a high contrast of 10 9 in polychromatic light (A/1 ^ 
0.02 /lo), corresponding to a 5 cr detection of 5 x 10" 9 . 

4. Conclusions 

We have numerically demonstrated that the SCC associated with 
a 32 x 32 DM enables us to detect Earths from space in a few 
hours when using realistic assumptions (zodiacal light, photon 
noise, read out noise, VLT pupil aberrations, and 20% band- 
width). SCC could be a good candidate to be implemented in 
the next generation of space telescopes. The technique involves 
two steps. We first use SCC to estimate wavefront errors and op- 
erate a DM that completes the correction in a few steps. To over- 
come the limitation linked to the DM uncorrectable high-order 
frequencies, we apply to the final iteration image, the SCC post- 
processing algorithm. This post-processing has yet to be opti- 
mized. 

SCC is one of the techniques under investigation for the E- 
ELT planet finder so-called EPICS. For this reason, we propose 
to consider the impact of different parameters, such as ampli- 
tude errors and turbulence residuals on the SCC performance. 
We will also test the compens ation for amplitude errors pro- 
posed by lBorde & Tr aub (2006). A preliminary study, which as- 
sumes a more realistic coronagraph (achromatic Four Quadrant 
Phase Mask), indicates that our algorithm converges but more 
slowly than with a perfect coronagraph. The quality of the ref- 
erence beam should not be important for SCC because of the 
filtering by the pinhole and the reduction in the beam diame- 
ter (D R ), which induces a wide diffraction pattern in the focal 
plane. Experimental validations of the SCC technique are also 
planned soon. 

We thank Pascal Borde and Anthony Boccaletti for useful 
discussions. 
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Under the same assumptions, we simulate four 2 x 10" 10 
companions at 1, 3, 5, and 1A /D (0.02, 0.06, 0.10 and 
0.14arcsec), including their photon noise. As shown in Fig.|U 
these Earth-like planets are detected in the SCC post-processed 
image after a total exposure time T of ~ 3 hours. The accuracy in 
the measured positions is a fraction of A/D (lower than A/ (2 D) 
). Fluxes are determined with a precision better than 20% for the 
three most off- axis companions. The coronagraph degrades the 
accuracy of the measured flux of the closest companion (lAo/D): 
the image of this companion appears to be slightly distorted be- 
cause the Earth contrast is just above the detection limit (Fig. [3} 
and a residual speckle is present at that position. The efficiency 
of the post-processing algorithm should be improved in future 
studies. We note that the correction area is larger in the fringe 
direction (from top-left to bottom-right) because of the residual 
chromatic dispersion effect (Eq.[2]). 

Similar results fo r high contrast imaging were demonstrated 
by iTrauger & Traubl (J2007I) in a laboratory experiment. They 



